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When the electromagnetic field is detected by stimulated emission, rather than by absorption, 
antinormally ordered photodetection can be realized. One of the distinct features of this photode- 
tection scheme is its sensitivity to zero-point fiuctuation due to the existence of the spontaneous 
emission. We have recently succeeded in experimentally demonstrating the antinormally ordered 
photodetection by exploiting nondegenerate stimulated parametric down-conversion process. To 
[ properly account for the experiment, the detection process needs to be treated with time-dependent 

. and continuous-mode operators because of the broadband nature of the parametric down-conversion 

' process and the wide spectrum of the pump that we used. Here, we theoretically analyze the antinor- 

I mally ordered intensity correlation of the continuous-mode fields by pursuing the detection process in 

, the Heisenberg picture. It is shown that the excess positive correlation due to zero-point fluctuation 

' reduces because of the frequency-distinguishability of the two emitted photon pairs. 

■ 

in : 

I. INTRODUCTION 

> ■ 

' In the context of the quantum theory of optical coherence Glauber found that the measurement operator for a 

■ photodetection can be associated with an annihilation operator Jj. His assertion rests on the fact that the detection 
<^ ' of a photon of energy huj is usually carried out with photoelectric effect, which absorbs the photon and dissipates 
C^l . energy hu; from the field. Any correlations measured by the absorption-based detectors are thus represented by the 
Op ■ expectation values of the normally ordered product of the annihilation and the creation operators Q. 

' The normally ordered photodetections can be said to be insensitive to zero-point fluctuations of the electromagnetic 

■ field because the photodetection probability for the vacuum state is zero, i.e., Tr[ |0)(0|a^a] = (0|a^d|0) — 0. This 
' insensitivity to zero-point fluctuations results in the photodetection process being logically irreversible^,^, that is, 

the premeasurement density matrix of the measured system not being uniquely determined by the postmeasurement 
-j-J ' density matrix and the outcome of the measurement. 

^ , When photons are, however, detected by stimulated emission, rather than by absorption, the detection responds 
^ ' not only to actual photons but also to zero-point fluctuations via spontaneous emission. Then, the antinormally 
ordered photodetection can be realized. Owing to the sensitivity to zero-point fluctuation the antinormally ordered 
r* 1 photodetection may provide an interesting alternative approach to continuously monitoring quantum system. This 
. kind of emission-based photodetector, what is called the quantum counter, was originally proposed by Bloembergen for 

detecting infrared photonsP|. The correlations measured by the quantum counters were then theoretically investigated 
by Mandel 5]. Ueda and Kitagawa showed that the quantum counter, unlike the standard photodetector, make it 
possible to perform a logically reversible measurement j^l . 

One way to realize the antinormally ordered photodetection is to exploit nondegenerate stimulated parametric down- 
conversion process. Using this process we have recently succeeded in experimentally implementing the antinormally 
ordered photodetection and measuring the antinormally ordered intensity correlations for coherent states @. The 
observed correlations were different from the normally ordered ones, as they showed excess positive correlations, i.e., 
the photon bunching effects, or the Hanbury-Brown— Twiss effectsQ- 

To properly account for the experiment, the detection process needs to be treated with time-dependent and 
continuous- mode operators Q because of the following reasons. First, in the experiment, a ultrashort pulsed field 
with the broadband spectrum was used as a pump of the parametric down-conversion to obtain large nonlinear re- 
sponse of the crystal and to overcome the slow response time of the photodetectors[30|. Second, as a nature of the 
parametric down-conversion, the spectrum of the down-converted fields themselves were broadband. Moreover, we 
used a spectral filter for increasing the coherence time of the down-converted field to be larger than that of the pump 
pulse in order to observe the field correlations. In the previous report 0, we only presented the simple single- mode 
treatment of the detection scheme. In this article, we give a theoretical treatment of the antinormally ordered pho- 
todetection process by using the continuous-mode field operators and pursuing the detection process in the Heisenberg 
picture. For the vacuum input, the resultant expression of antinormally ordered correlation is equivalent to that of Ou 
et aL|lli|. which is obtained by treating the continuous-mode parametric interaction in the Schrodinger picture. Our 
expression for more general coherent field input reasonably explain our experimental result reported previously (6J. 
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II. SINGLE-MODE TREATMENT 

The original theoretical proposal of the antinormally ordered photodetection, i.e., the quantum counting^ as well 
as the subsequent theoretical analyses^ Q utilized the cascade emission process of atomic system. This resonant 
configuration based on the atomic cascade emission can be translated into the non-resonant configuration based on the 
parametric down-conversion. Our experimental implementation of the antinormally ordered photodetection followed 
this line. However, to measure the antinormally ordered correlation is by no means a straightforward task because 
all normally ordered correlation terms should be eradicated. Now, we begin by showing how our devised methodlsl 
works for measuring the antinormally ordered intensity correlation. Here, we will treat the parametric down-conversion 
process in the simplest way and thus treat only a few single-mode operators relevant in this measurement process. 
In section fHIl we will analyze the measurement process with the time-dependent, continuous-mode operators. A 
schematic illustration of the antinormally ordered intensity correlator is shown in Fig. ^ Now, let us track the time 
evolution of an annihilation operator, ain, which represents the signal field. First, operator aln is coupled with an 
operator, bin, via the parametric interaction with a monochromatic pump field. The coupled spatial differential 
equations for a mode of the signal field, a(z,a;), and a mode of the auxiliary field, b(z,LUp — lu), in the nonlinear 
dielectric medium with the monochromatic pump field of frequency LUp are written as 



(1) 



= -s(LUr, - uj, uJ e-'''('^p-'.-'p) b(z, UJ„ - LU) g-"^''^'^^ -'^^^ . (2) 

dz 

where the phase-matching parameter, Ak{ujp,uj), is defined by 

Ak[cjp,uj)= — . (3) 



Here, , rt^, and -cj are the refractive indices of the nonlinear dielectric medium at pump field frequency ujp, 
signal field frequency to, and auxiliary field frequency ujp — to, respectively. s{uj,ujp) and d{Lj,ujp) are explicitly given 

byiS 



and considered as the coupling constant between the two modes, a{z,uj) and b(z,ujp — lu). Here, F(ujp) and 4>{ijJp) 
correspond to the power spectral density and the phase of the pump field, respectively, x*^^' (tj, tjp) is the second-order 
nonlinear susceptibility. We can easily show that s{uj,ujp) e^^'^'^^'^p^ = s{u!p — uj,ujp) e^'^^'^''~^'^''\ The two modes, 
a{z, uj) and 6(z, ujp — u), as well as the pump field are assumed to be propagating almost colinearly in the z direction, 
here. When the phase-matching condition is perfectly satisfied, i.e., Ak{ujp,uj) — 0, we have the following solutions: 

Qout = a{L,oj) ~ Oin cosh [sL] — bin e''^ sinh[sL] 

= sUOa^nSabiO (5) 
bout = b{L, ujp — uj) — bin cosh[sL] — aln^e^"^ sinh[sL] 

= SliOblnSabiO (6) 

with the initial conditions a{0,uj) = aln and 6(0, Wp — w) = bin- Here, L is the length of the crystal. Sab is the 
two-mode squeezing operator, which is given by 

SabiC) = explCalnbin - Ca'iJbin] (7) 

with C = sLexp[it9]. Intuitively, one photon in the pump field is split into two photons in such a way that the total 
energy of the fields is conserved during the process. An energy diagram of this process is shown in the inset of Fig. ^ 
Next, by dividing the field represented by bout into fields bi and &2 with a half-wave plate and a polarizing beam- 
splitter, the standard Hanbury-Brown— Twiss interferometer |7ij for field bout is formed as shown in Fig.^ Here, the 
imperfect detectors and several optical losses are taken into account by introducing auxiliary beam-splitters with 
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FIG. 1: Schematic illustration of the antinormally ordered Hanbury-Brown— Twiss-type interferometer based on stimulated 
parametric down- conversion. The inset shows an energy diagram of parametric process. 



vacuum fields vi and V2^. Output fields di and of Fig. are then written as 

di = .Jrfl [Tbout + Tlv] + i^l - 771 vi, 

d2 = y/m [R-bout + Tv]+i^l-ri2 W2, (8) 

respectively. Here, 771 and 772 are the total photodetection efficiencies for detectors 1 and 2, respectively; TZ and 
T denote the reflection and transmission coefficients at the beam-splitter, respectively, and can be varied with the 
half-wave plate and the polarizing beam-splitter; and v is an auxiliary vacuum field introduced from the empty port 
of the beam-splitter. 

Since the modes relevant to operators &i„, v, ui, and V2 are initially vacua, the surviving contribution to the 
number of photodetection events of the normally ordered photodetection for fields di and d2 coincides with that of 
the antinormally ordered photodetection for field aln up to a constant factor, i.e., 

(ndi) = {di di) ^ rii\T\^ s\VLh^[sL\{a'inaiJ), 

{fid^) = (dVrf2) = f?2|7?.psinh^[sL](ai„ai„'''), (9) 

respectively. Furthermore, the number of coincidental photodetection events of fields di and ^2 (i-e., the interbeam 
correlation for the two fields) results in 

{fid^fid^} = {di did2 d2) = ??i?72|'?'n7?.|^ sinh''[sL](a;„a;„aJ„''' (10) 

where we use the commutation relation for each operator and relations TTZ* = —TZT*, which should be satisfied 
when the lossless beam splitter is used. Thus, the surviving contribution to the coincidences turns out to be made 
only by operators a^n and ai„^ in antinormal order. Consequently, we can evaluate the antinormally ordered intensity 
correlation for field din as follows: 

(2) _ {fldi'hd2) _ {O'inO'inO'in^ 0,171^) _ {2[A]) /i -i \ 

\ndi)\nd2) {ainttin^ ) [amain ) 

Note that Eq. holds regardless of the splitting ratio at the beam-splitter, the quantum efficiencies of the detectors, 
and the optical losses. In this way, we can realize the measurement of the antinormally ordered intensity correlation 
by exploiting the parametric interaction, Eqs. © and 

III. CONTINUOUS-MODE TREATMENT 

To analyze the correlation observed in our experiment, the detection process should be treated with time-dependent 
and continuous-mode operators Hj. In this treatment, the measured antinormally ordered intensity correlations in 
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Eq. becomes time dependent. Moreover, since the response-time jitter of the detector is larger than the pump pulse 
duration (the duration of the parametric interaction) but smaller than the time interval between two successive pulses, 
the relevant information on the time dependence is embodied in the integrated number of the delayed coincidences 
over the response-time jitter of the detectors 9, 10]. In the following we treat the antinormally ordered photodetection 
process by using the continuous-mode field operators and pursuing the detection process in the Heisenberg picture. 
For the vacuum input, the resultant expression of antinormally ordered correlation is equivalent to that of Ou et 
a/.|ll|. which is derived by the continuous-mode analysis of the parametric interaction in the Schrodinger picture. 

We start by reconsidering the quantization procedure to deal with the realistic photodetection process, in which 
a photon in a travelling field with a broadband spectrum is dissipated at a detector's photosensitive surface (cross- 
section A), during period T. Then, we tackle the continuous-mode parametric interaction in the Heisenberg picture. 
Finally, the antinormally ordered intensity correlations of continuous-mode fields are derived. 



A. Quantization procedure for travelling field 

Following the method developed by Huttner, Serulnik, and Ben-Arveh[l2|. we review a quantization procedure 
for treating the travelling field with multimode operators. An annihilation (a creation) operator, which annihilates 
(creates) a photon of frequency LOm travelling through plane z — zq with cross-section A during period T, is defined as 
a{zQ, ojm) (a^(-zo, t^m))- Here cum = 2'Km/T . The positive-frequency parts of electric and magnetic fields in vacuum are 



written as E^+\zo,t) = Y.m \j ^^(^izo.oo^)e~'^"^' and B(+)(zo,0 = Y.m \j 2^^a{zo,u}^)e-''^^\ respectively, 

where Ti is the Planck's constant, eo is the electric permittivity of vacuum, and c is the speed of light in vacuum. Then 
the normally ordered Poynting vectorfsLIT^ takes the form 



^(zo,t) = c2eo[^(-)(zo,i)S(+)(zo,<) + ^(-)(zo,t)^W(zo,t)] 

n 

2AT 



^Tr V^™^"' [a\zo,ujm)a{zQ, )e't"'" ""'1* -f ff.c.]. (12) 



Here and throughout this subsection, we exclusively use () for denoting variables to avoid any confusion. By integrating 
the Poynting vector over the period, T, we have the total energy. 



A I dtS{zo,t) ='y^hLUrna,\ 
■Jo 



Zo,U!rn)d{zo,U!rn), (13) 

which flows through the plane, z = zq, with a cross-section, A. The quantization volume for the operators a(zo,Wm) 
and {zq,uj„i) is thus AcT. It is important to note that the total energy, Eq. H13|l . which is calculated without 
assuming a real cavity, does not involve zero-point energy. Thus the total energy, Eq. (|13|l . can be considered as the 
energy that is dissipated at a detector's photosensitive surface at z = zq with cross-section A, during the period, T. 
Here, the equal-space commutation relations (ESCR): 

[d{zo,UJjn),a^{zo,UJjn')] = Sjn,m' (14) 

are introduced T^]. It can be verified that these ESCR are consistent with the usual equal-time commutation relations 
(ETCR) of the canonical quantization procedure and are even valid for the case of a field in a dielectric medium [12i . il3l| . 
The number operator of photons localized in the finite volume, AcT, can be then written as 

N{zo) ^y2d^{zo,t^m)d{zo,ujrn) ^ / dti{zo,t), (15) 

where the operator, I{zQ,t), represents the photon flux at position z = zq and cross-section A. The photon flux 
operator, I{zQ,t), takes the simple form: 

iizo,t) ^Sizo,t)d{zo,t). (16) 

with the photon flux amplitude operator, 

d(zo,t) = -j=a(zo,uj„,)e-''^"^'. (17) 
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Here, we stress that even though the number operator, Eq. (|15|l . is taken at position z = zq, it is not the number of 
strictly locaUzed photons a,t z = Zq, rather it means the number of delocahzed photons over the length, cT, and the 
cross-section. A; otherwise a localized number operator for photons would contradict the well-known fact that there 
does not exist a probability density for the position of photons [T^ Hsl . 

The continuous- mode extension of the above multimode treatment is straightforward by taking the limit, T — > oo, 
with almost-always-valid narrow-bandwidth approximation, by which the range of integration over frequency lu is 
extended to cover from — oo to oo 8J. For instance, the photon flux amplitude operator, Eq. H17|l . becomes d{zQ,t) = 

Hereafter, we will work with the continuous-mode extension of the annihilation and the creation operators, which 
satisfy the continuous-mode ESCR, 

[a{zo,Lo),d^zo,u;')]^S{Lu,Lu'), (18) 

which is analogous to the multimode ESCR, Eq. H14|) . To include the effects of the spectral filter with the amplitude 
transmission function, ^(w), we can further modify the photon flux amplitude operator, Eq. (|17|l . as 

1 

d^{zo,t) = ^ dLoi{uj)d{zo,uj)e-'^K (19) 



/27r 

Then, the number operator, Eq. (|15|l . is changed to 



N{zo) = I dti{zQ,t)^ / dtd^zo,t)d{zo,t) 

duj\i{uj)\' d\za,uj)d{zo,uj). (20) 



The operator, N{zq), can be identified as a positive operator-valued measure\\J\ for a photodetection measurement 
with the spectral filter characterized by ^(w) during infinite time. This expression is an approximation of the rather 
slow response time of the standard photodetector and indicates that the intermode coherences disappear during the 
photodetection0| . 



B. Continuous-mode analysis of parametric down-conversion in the Heisenberg picture 

In Sec. ^ the analysis of parametric interaction was carried out under the perfect phase-matching condition, i.e., 
Afc(wp,a') = 0. Besides, the pump field was taken as to be monochromatic. In our experiment 6], however, the 
condition, Afc(wp,w) = 0, did not perfectly satisfy in generating the down-converted fields, and the pump field was 
a ultrashort pulsed field with the broadband spectrum. In this subsection, to appropriately examine the experiment 
we will develop a theory of the stimulated parametric down-conversion process for the continuous-mode fields in the 
Heisenberg picture. 

To elaborate the theory of the parametric process so as to include the pump field spectrum and the frequency spread 
of the down-converted fields, let us go back to the two-mode squeezing operator, Sab of Eq. 10 . The continuous-mode 
expansion of the squeezing operator can be given by 

^, = exp[P;-P;], (21) 
where may be called the frequency- entangled photon-pair creation operator^ which is defined by 

/oo 
dwidw2l{i^l,Ul2)(i)(lOl,UJ2)alJ {uJi)brn {^2), (22) 
-00 

where 

C(wi,W2) = s{ijJi,bJi + W2)-^^exp[ii9(wi,a;i -1-^2)]- (23) 

Since the refractive indices, n^j^+i^^, n^j^, and rij^^, the frequency multiple, (wi -f 1^2)^1 ^^2, and the second-order 
nonlinear susceptibility, x*^^' (a;, cjp), in s{ijJi, loi 4- 0^2), which was given by Eq. can be approximated by constants 
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with respect to frequencies wi and oji over the spectral filter's bandwidth (which is determined by the profile, 
in Eq. (|20|l '). they can thus be factored out as 



C(wi,w2) « Cy^pK + ^2). (24) 

The remaining frequency-dependent terms of Pj in Eq. (|22|l are the pump field profile, Fp{uJi + UJ2), and the phase- 
mismatching function, (/)(Li;i, 072)- Then, the continuous-mode counterparts of the down-converted fields, Eqs. © and 
®, can be written as 

aout((^) = Sl{()ain{uj)Sc{C), (25) 
Ct(^) = 5l(C)C(c^)^c(C), (26) 

respectively. Here, SdC) is defined by 

Sc{C)^exp[CP^-CPl], (27) 

with 

P| = / dwi dtj2 i/Fp (tJi + W2 ) (/>(wi , ^2 ) aln H^i) hn {^2 ) 

da;idcj2<&(wi,t^2)ain^(wi)6m\w2)- (28) 



C. Antinormally ordered intensity correlation 

By incorporating the results of two subsections above, we will now show how the antinormally ordered intensity 
correlation of single-mode field given by Eq. pf |l is altered when the field is replaced with continuous-mode field. 
First of all, the output operators, di and 6,2, of Eq. (|HJ are rewritten as the flux amplitudes, 

diit) = —= / duj(,{uj)di{Lo)e-"^\ 



1 



d2{t) = ^ / da;C(w)d2(^)e"""*, (29) 



with 



M^) = Vm [nbout{0j)+T v{Lu)]+i^l-Tj2V2{i0), (30) 

respectively, where bouti^) was given by Eq. (|26|l . Here, the parameters, T, 7?,, 771, and 772, are assumed to be 
independent of the frequencies of the fields. As noted in Sec. ^ the modes relevant to operators v{uj), wi(w), and 
V2{lij) are supposed to be initially vacua. Then, the number of the photodetection events for the field di{uj) with 
detector 1 reduces to 

{Nd,)^ dt{di {t)Mt)) = — dt duodLu'C{io)i{oj'){di{Lu)ddLu')) 



In 

1 f°0 j'OO ^ 

= — J7i|T|2/ dudu'C{^)a^'){b:j{Lo)b:ut{io'))e'^^-^^' 
J —00 J —00 

1 /-oo j-OO 

= —in\T? dt dc^dc.'r(c^)^(c^')(^l(C)fe.n ('^)5c(C)^KC)fc.«(^')5c(C))e^^"-"^* 

J -oa J -00 

= 2^^il^l'/ dt j dwdJC{oj)a^%Sl{()bJ{oj)bUu')S,{0)e'^^-'-^\ (3f) 



00 —00 



where we use Eq. (|26ll . Furthermore, since the modes relevant to bin are also supposed to be vacua |0);,, we have the 
following expression for &i„(w')S'c(C)|0)f, to the first order in (: 

bU^')s,{C)\o)b - bUuj'){i + [CP* - CPI] + liCP-P - cpI? + ■ ■■)\o)b 

« -CbU^')Pl\0)b 

f°° „ f 

f°° , f 

^ — OO 
POG 

= -C / dcu[<i>{uj[,ij')a2j{cu[)\0)b. (32) 



Then, Eq. ^ leads to 



OO poo 



(Na,) - ^ml^nCI' / dt I dcjdij'du;idu;[C{^)a^') 



— OO 'J —OO 



x$*(c^i,c^)<i>K,c^')(«m(^i)a;„^(wi))e*("-'^')* 
= m|rp|C|2 / diodLOidLo[\aiu)\^^*iu;uiuMu;[,u){a2„iLOi)a2n\L^[))- (33) 



Likewise, the number of the photodetection events for the field d2{LLj) with detector 2 can be given by 

poo 

(iVd,) =r,2|7^p|Cp / dcc;da;ida;i|e(a;)|'$*(a;i,w)$(w;,c^)(a;„(a;i)a-„^K)). (34) 



Therefore, the normally ordered photodetection for field di{uj) and that for d2{uj) again coincide with that of the 
antinormally ordered photodetection for the input signal field ai„(a;), although, this time, the frequency-dependent 
factor, ^^{uji,uj)^{uj[,uj), is added. 

The number of the coincidental photodetection events in detector 1 and detector 2 is given by 

{Nd,,d,) = / dtdt\di {t)d^{t)d2 {t')d2{t')) 



-| pOO pOO 

= — dtdt' dujduj'duj"duj"'C{i^K{^')Cioj")£.{uj"') 

277 

x(dV(cj)rfi(a;')dV(t^'')d2(tc''''))e'["""'l*+'t"''""'''i*'. (35) 

With the vacuum conditions of the modes, v{uj), vii^^), and V2{lo), the surviving terms in the expansion of Eq. H35|) 
are found to be 

1 pOO POO 

{Nd„d,) = ^viV2 / dtdt' / di,duj'duj"duj"'c{i^n^')c{i^")a^"') 



— oo J —oo 



27r 

xe^["-"'l*+'["''-"'''l*'(|r^7^|2(C/(u;)5;„i(t^0C^\c.'0Ct(^''')) 

1 poc poo 

= —ViV2 / dtdt' / dujdLu'duj"duj"'C{(^Mi^')C{i^"M^"') 

^TI" J — oo J —oo 

xe^["-"'l*+'l"''-"'''l*Vn7^p((6;„J(t^)5;„t(t^0CJ(cc''0Ct(c^''')) 

-{bout {uj)bout{'^"'))5{uj' - w")) 

1 poo poo 

= i^Tlim\T'?\T^? dtdt' dudu'dLu"doj"'C{^)^{i^')C{^"W') 



27r 

x(6;.J('^)CJ(c^")Ct(a;')Ct(c^"'))e'^'"""'^*+''"""""'^*' 
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-| p oo poo 

= T^?7l?72|T|2|7^|2 / dtdt' dudu'dLu"doj"'C{^)^{i^')C{^")^W") 

J-ryO J-OO 

x{Sl{Obl\io)b:J{co'%„{Lo')b:r.ico'nSc{0)e'^^^^^ (36) 

where we use the relation, TTZ* = —T*TZ. Since the modes relevant to bin are also supposed to be vacua |0)fc, we 
have the following expression for bin{iLi')bin(oj"')Sc{O\0)b to the second order in C: 

bULj')bUu;"')Sc{Cmb 

= + [CP^ - CPl] + l[CP<!> - (PI? + ■ ■■)\o)b 

« ^CHUuj')bUco"'){Plmb 

-(^ b^n{uj')bin{oj"') / duj[duj'2duj'{'duj!^'<S>{uj[,Uj'2)<i>{uj'l',Uj'^') 



2 

xa'iJ{uj[)b,J{u2) all ^ ( w"' ) ^ Wj" ) 1 0) 6 

1 

^ J— oo 

X {d{uj' - u;!,)S{uj"' - J^') + - ij'^')S{uj"' - lo'^)) \0)b 

\e r doj'.doj'i'mLoi^'MLo'i',^'") 



where we use the relation: 



Then, Eq. ^ leads to 



+<i>K,c."')<i>K',c.'))ai„^K)a;„nO|0)6, (37) 

[bUio')bUuj"'), bj{uj'^)bj{uj'^') ] |0)6 
= {6{lu' - lo'2)S{lu"' - uj'^') + S{lu' ~ lu'^')S{lu"' - to'^)) \0)b- (38) 



{Nd,d2) = ^^r^l772|r^7^|2Al:lJ-^ / dtdt' dujduj'duj"duj"' dujiduj'.duj'ldcu'/' 
(27r)^ 4 



' —OO 



(2^)2 

xr(c.)e(t^')r(^")c(c^"') 

X ($*(wi, t^)$*(tji', Oj") + $*(cJi, tj")$*(Wi , t^)) 

x(a;„(a>i)a-„(c.;')ajJ(c.;)a-„t(c.'/'))e'^['^-'^'l*+^[""-""'l*' 
= ^7l?72|r^7^|2^^L^ / dwdto"dtoidto[dto'{dto'{'\^{iu)\^\^iu;")\'' 

X ($* (wi, t^)$* (wi', t^") + $*(tJi , cj")$*(t^i , t^)) 
X {^{lu[,u)^{lu';', lo") + $(cc;i, c<j")$(c^i', w)) 

X (ai„(a;i)a-„(a;")ai„^(w'i)ai„'''(cj"')). (39) 

As in the case of Eq. fTIII , the surviving contribution to the normally ordered coincidental photodetections for fields 
di(w) and d2{Lo) turns out to be made only by the input signal field alnito) in antinormal order, though the complicated 
frequency-dependent factor is added in this continuous- mode counterpart. An attention to the frequency response of 
the detection (parametric interaction) needs to be paid when the polychromatic-light detection is to be concerned. It 
should be emphasized that this sort of situation also occurs in measuring the normally ordered correlations jl9j. 

Finally, we apply the above continuous- mode treatment to our experiment 0, in which the antinormally ordered 
intensity correlations for coherent states were measured by utilizing stimulated parametric down-conversion. Let the 
input signal field be a continuous-mode coherent state|^] with spectrum amplitude a{uj): 



\a)a = exp[-- / dw|a(w)|2] ^ ^-^^ ^^|0)a. (40) 

n=0 
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By rewriting Eqs. (|33|l . (|34() . and H39|l in the normal order, we can exploit the relation, ain{Lo)\a)a = oi{^)\o)a- The 
resultant expressions are respectively written as 



J —OO 

x<^*{uji,uj)<^{ijj[,uj){a*{uj[)a{uji) + 5{uj[ - uji) 

/OO 
dL0dL0ML0)\^\^{L0l,iu)Wn+l) (41) 
-OO 

/OO 
-OO 

x<^*{uji,uj)<^{uj[,uj){a*(uj[)a{uji) + 5{uj[ - uji) 
= 7,2|7^^CP / dcjdcJi\auj)\^Mu;uUj)\^n+l), (42) 



and 



(iVdi.d.) =^7l^2|T|2|7^|2^^ / da;dc."dc.idc.idc.i'd<|C(a;)ne(c^")l' 

^ ^ — OO 

X ($*(tJi, t^)$*(tji', u") + $*(tJl, cj")$*(t^i , t^)) 
X {^{lj[,uj)<^>{uj'I', uj") + $(tj;, cj")$(t^i'. t^)) 
X (a*{uj'l)a*{Lu';')a{oji)a{Lu[) + <5(cji - Lo'l')a* iuj'l)a{oj[) 

+6{uj[ - tji")a*(u;i)aK) + H^i - )a*K')a('^i) + '^(t^i - - w"') 



/OC 
da;da;"dc.idc.i'|f(c.)n^(^")P 
-OO 

X ( I $ (a;i , c.) n $ (a;i' , c." ) P (n' + 2n + 1 )) 

+$*(wi, w)$*(wi', w")«'(wi, w")*(wi , w) (2n -I- 1))) . (43) 

Here, for simplicity, we approximate a* (LO\)a(L^'^ by n6(uj\ — wj.)? where n describes the average photon number of 
the input signal field, and likewise the other pairs such as a* {Lo'i)a{Lo'{") and a*(w")'^('-^i')- These approximations 
correspond to regarding the input signal field a\n{'^\) as a stationary and frequency- independent field and may be 
justified by the fact that the coherence time of the input signal is much longer than that of the down-converted 
fields and the frequency-dependent characteristics of the input field can be pushed into that of the pump field, i.e., 
Fp{LO\ +^^2)- Then, we have the modified antinormally ordered intensity correlation for the continuous- mode field, 

■ 



where 

\^Ji^dJ'do,^dJi\i{uj)f\i{JT^*{u:^,u:W 

7= 5 • (45) 

(Oc.da;i|C(^)P|<J>(c.i,..)P)^ 

Here, 1 -t- 7 is exactly equivalent to the form of intensity correlation g^^^ ^ for the spontaneous parametric down- 
conversion, which was derived by Ou et aL|ll|. Physically, value 7 can be viewed as a measure of the indistinguisha- 
bility of two emitted photons responsible for a coincidental photodetection. By using a narrow spectral filter, we 
can eliminate the frequency correlation of the down-converted fields, which gives rise to the distinguishability of two 
down-converted photon pairs and the reduction of value 7. Then, since the functions, ^(lo\^uj) etc., can be factorized, 
the correlation, Eq. H44|) . reduces to Eq. _ 

Now, let us move on to the analysis more specific to our experiment. The experimental results Q were indeed in good 
agreement with the form, Eq. 14411 with 7 = 0.45. To check whether this value of 7 is really reasonable, we will follow 
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the recent analysis due to de Riedmatten et a/.|20|, in which the multiple spontaneous parametric down-conversion 
process was treated. First, for ease of calculation, we approximate the filter's spectral profile, by 

and the function, $(aji,W2), by 

Hco^u:,) « ^F.ico.+u.,)^ -4-^exp[- ^"^+;^7"^)' ], (47) 

where A|, and are the variances of the power spectral density for the filter and the pump field, respectively, 
and and Op correspond to the center frequency of the filter and pump, respectively. Here, phase-mismatching 
function, 0^2)7 is assumed to be one since the frequency spread of the down-converted fields might be essentially 

determined by the filter's spectral profile, in our experiment. Since both functions, |<^(a;)p and [^(ijJi, Li;2)P, 

were assumed to be the normalized Gaussians in Eqs. H46|l and H47f) . the denominator of Eq. (|45|l becomes one after 
the integration. The calculation of the numerator in Eq. H45|) is, on the other hand, rather tedious. Using square 
completion technique, we can integrate (f)*{uji,Lu)(f){u;i,U!") with respect to uji and (/)* (loI , Lu")(f){uji , u) with respect to 
w"- Then, we have the following reduced form of 7: 

1 up' "'^ 1 " 

^ " 2^ J '^'^^'^"°'^P["^]^^P[^]^^P[-iA2^^' " (48) 

where, without loss of correctness, we set i7 = 0. Using square completion technique again, we first integrate Eq. 148(1 
with respect to u" and then with respect to uj. Then, we get 



27rA2, 



\/27rr2\/27rr'2, (49) 



where j4- = ^ + and ~ + Thus, the final simplified form of 7 can be written aspof 

7= (50) 

In our experiment, we used the 5-nm-FWHM filter and the 10-nm-FWHM pulse for the pump field. Thus, we have 
7 « 0.9 from Eq. jSUJ). Besides, since the the filter's spectral profile was found to be not Gaussian, rather Lorentzian, 
thus the coherence time of the down-converted fields stretched by the filter might be less than that achieved by a filter 
with the Gaussian spectral profile as was assumed to be in Eq. 1(501) . In addition, the pump pulse duration might be 
broadened in the course of propagation. Taking account of these factors, value 7 further reduces by about 10%. The 
value is still higher than the observed value of 0.45. However, Eq. ((50(1 was derived by concerning only the temporal 
mode mismatching, i.e., the frequency dependence of the down-converted fields. The imperfect alignment, i.e., the 
spatial mode missmatching, may give rise to the further reduction of the value, 7. In fact, the coupling efficiency 
between the down-converted field and the single-mode fiber was 60% on average. Roughly speaking, this imperfection 
reduces the value, 7, by 60%^. Thus, the experimental value, 7 — 0.45, may be the reasonable end result. 



IV. CONCLUSION AND OUTLOOK 



We have analyzed a scheme of antinormally ordered photodetector, i.e, quantum counter, which utilized a stimulated 
parametric down-conversion process, with time-dependent and continuous- mode operators. It has shown that the 
operator ordering of the intensity correlations measured by such detectors is antinormal even when the measured field 
has the finite spectrum, though the correlation in general include frequency-dependent factors. The experimentally 
measured antinormally ordered correlations can be well explained by the present analysis. 

The emission-based antinormally ordered photodetection may provide an interesting alternative for monitoring 
quantum systems owing to its sensitivity to zero-point fluctuations. The noise reduction in the optical signal 
amplification 21], and the realization of the logically as well as physically^, |2^ l23j | reversible measurements may 
be just a few examples of the possible applications of the antinormally ordered photodetection. 
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